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ABSTRACT 

An integral part of the Unified Model for Active Galactic Nuclei (AGNs) is an axisymmetric 
obscuring medium, which is commonly depicted as a torus of gas and dust surrounding the central 
engine. However, a robust, dynamical model of the torus is required in order to understand the 
fundamental physics of AGNs and interpret their observational signatures. Here we explore self- 
similar, dusty disk-winds, driven by both magnetocentrifugal forces and radiation pressure, as an 
explanation for the torus. Using these models, we make predictions of AGN infrared (IR) spectral 
energy distributions (SEDs) from 2 — 100 /im by varying parameters such as: the viewing angle 
(from j = 0° — 90°); the base column density of the wind (from Nnfi — 10^^ — 10^^ cm~^); the 
Eddington ratio (from L/LEdd — 0.01 — 0.1); the black hole mass (from Mbh — 10^ — 10^ M©); 
and the amount of power in the input spectrum emitted in the X-ray relative to that emitted in the 
UV/optical (from Uox = 1.1 - 2.1). We find that models with Nh,o = 10^^ cm^^ L/LEdd = 0.1, 
and Mbh > 10^ Mq are able to adequately approximate the general shape and amount of power 
expected in the IR as observed in a composite of optically luminous Sloan Digital Sky Survey 
(SDSS) quasars. The effect of varying the relative power coming out in X-rays relative to the 
UV is a change in the emission below ~ 5 /im from the hottest dust grains; this arises from 
the differing contributions to heating and acceleration of UV and X-ray photons. We see mass 
outflows ranging from ~ 1-4 M© yr~^, terminal velocities ranging from ~ 1900-8000 km sec^^, 
and kinetic luminosities ranging from ~ 1 x 10^^-8 x 10''^ erg s~^. Further development of 
this model holds promise for using specific features of observed IR spectra in AGNs to infer 
fundamental physical parameters of the systems. 

Subject headtngs: galaxies: active — galaxies: Seyfert — hydrodynamics — MHD — infrared: general 
— quasars: general — radiative transfer 
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1. Introduction 

Active galactic nuclei (AGNs) are often set 
apart from other astrophysical objects by their 
powerful spectral energy distributions (SEDs) that 
radiate a substantial amount o f power across th e 
electromagnetic spectrum (e.g., Elvis et al]ll994 ). 
In addition to being powerful, AGNs are also a 
motley lot, known for displaying a wide range of 
observable characteristics among them. Objects 
of two types - radio-loud and radio-quiet - have a 
range of intrinsic luminosities, and may display ei- 
ther only narrow, o r both narr ow and broad spec- 
tral emission lines (| Antonuccil 1 1 993i) . 

1993t 



The Unified Mode l (e.g., lAntonucci 



Urrv fc Padovanilll995l ) is a phenomenological de- 
piction of AGNs that posits that the nature of 
many of these differences, particularly the pres- 
ence or absence of broad emission lines, can be 
explained primarily as a function of viewing an- 
gle. Pictured at the core of this model is the 
"central engine," composed of the supermassive 
black hole (with a mass from ^ 10^-10^ Mq) and 
its sub-parsec-scale accretion disk, where matter is 
heated as it falls towards the black hole. Another 
key feature of the Unified Model is a structure 
that surrounds the accretion disk: an axisymmet- 
ric, obscuring medium, which covers many lines of 
sight to the accretion disk around the central black 
hole. This obscuring medium is commonly de- 
picted as a torus of optically t hick dusty gas which 
can e xtend up to 100 pc fe.g.. lAntonucci fc Miller 
19851 ). This dusty torus, which is also the main 
source of infrared emission in an AGN, can block 
the accretion disk continuum and broad emis- 
sion lines produced within the Broad Line Region 
(BLR) by obscuring them with du sty gas (e.g., 
Sanders et al]ll989l : lBarvainislll990[ ). 

A consequence of this geometry is that when the 
AGN is viewed face-on, we observe a "Type 1" ac- 
tive galaxy displaying both narrow lines (permit- 
ted and forbidden) with line widths of hundreds 
of km s~^, as well as the broad lines emanating 
from the exposed hot (~ 10^ K) and bright pho- 
toionized gas around the central accr etion disk 
with widths of thousands of km s~^ ( Antonucci 



19931 ). When viewed edge-on, the broad lines are 
obscured by the optically thick dusty torus, and 
we observe only the narrow lines generated on 
much larger scales characteristic of "Type 2" ac- 



tive galaxies. 

As a further clue to the structure of the torus, 
there is observational evidence for a luminosity de- 
pendence of the fraction of Type 1 and Type 2 
AGNs. For exampl e, in a deep hard X-ray survey, 
Steffen et akl (|2003l ) show that Type 1 AGNs dom- 
inate at higher luminosities, while Type 2 AGNs 
become an important component of the X-ray pop- 
ulation o nly at low e r, Sey fert-like luminosities. 
Similarly, iHao et al.l (|2005al ) shows that the frac- 
tion of Type 1 AGNs increases with lum inosity of 
the [ OIIII narrow emission line (see also ISimpson 
20051 ). and studies using radio-selected AGN sam- 
ples consistent l y pq ii it to the same trend (see 



Hill et all Il996t ISimpson fc Rawhngsl [2000t 
20041) . 



Grimes et al 



Assuming the inner wall of the torus is set at the 



radiu s at which dust sublimates (e.g.. lElvis et al 
1994I ). then that distance is farther out in lu- 



minous quasars (the more luminous relatives of 
Seyferts). If the height of the torus for AGNs of 
varying luminosities remains approximately con- 
stant, then the opening angle of the torus must in- 
crease with luminosity. This leads to an expected 
dependence on luminosity in the observed f r action 
of Type 1 and Type 2 AGNs. iLawrencd (|l99l[ ) 



first suggested this idea, known as the "receding 
torus" model (though see iLawrence fc Elvisll2oTol 
for a different perspective). It is not clear why 
the height of the torus would be constant in all 
objects, but nevertheless a luminosity dependence 
is expected, as it would be highly unlikely for the 
height to vary in suc h a way as to p roduce a con- 
stant opening angle ( SimpsoiJ 20051 ) . 

This highlights the need for a more fundamen- 
tal understanding of the torus. As the structure 
at the interface between the AGN and its host 
galaxy, the torus plays a vital role in the make-up 
of a quasar. Yet in the past, it has typically been 
mode l ed as a stati c struc t ure (e.g., Pier fc KrolikI 



1992; Fritz et al. 2006: Nenkova et al. 



2008 



Schartmann et al.i2008. : iHonig fc Kishimotoll201 



which is likely only a rough approximation to the 
active inner region surrounding an accreting black 
hole. Indeed, for a geometrically thick torus to 
be supported hydrostatically, it must have a tem- 
perature on the order of 10^ K or higher; dust in 
such a hot, high-densi ty environment would most 
certainly be destroyed ( Dullemond fc van Bemme] 
However, dust signatures such as the "in- 
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frared hump" are observed ( Sanders et al. 1989I) . 
along with siUcate emission and absorption fea 



tures (e.g., ISiebenmorgen et al 
2005bl) . Therefore, dust near 



2005: Hao et al 



the sublimation 
radius is certainly present in a wide variety of 
AGNs. One way of getting around this inconsis- 
tency is to generate a static, geometrically thick 
torus through infrared radiation pressure (e.g.. 



Krolild 120071 : IShi fc KrolikI l2008f ): however, this 



treatment neglects the effect of radiation pressure 
from UV photons. 

One way to help us understand the structure 
and physics of the dusty torus is to examine the 
IR regime of quasar SEDs. Dust reprocesses a sig- 
nificant fraction of the accretion power in the UV 
and X-ray, and re-emits it in the infrared, acting 
as a sort of bolometer for the quasar. Indeed, the 
1 - 100 fitn IR "hump" of a quasar spectrum ac- 
count sjornea£ly40%of_tlie^ luminosity 
fe.g.. lElvis et al.lll994 [Richards et alll2006[ ). As 
commonly seen with Spitzer, AGN IR spectra also 
show specific structures, such as 10 ^m silicate 
features, that any successful model must explain. 
While Type 2 objects commonly show silicate in 
absorption, emission is often observed in Type 1 



object s (e.g. ISiebenmorgen et al.|[2005t Hao et al 
2005bl) 



High spatial-resolution imaging adds further re- 
quirements to any successful torus model. Using 
the Very Large T elescope Interferometer (VLTI), 
Jaffe et al.l (j2004l ) obtained interferometric obser- 
vations of NGC 1068 from 8.0 to 13.0 ^lm that 
show different dust temperatures found within 
close p roximity; this hints at st ructure within the 



torus (jSchartmann et al.l 120051 ) . Similarly, more 



recent VLTI observati ons (from 8.0 to 13.0 u m) of 
the Circinus AGN by iTristram etall (|2007l ) show 
irregular behavior (such as fluxes with rapid an- 
gular variations) that can be explained by dumpi- 
ness. Further, they find that the radial tempera- 
ture dependence of the dust indicates that some 
of the outer dust is exposed directly to radiation 
from the nucleus. Together, these observations in- 
dicate that the dust distribution in the torus is of 
clumpy or filamentary composition. 

Previous researchers have tried different ap- 
proaches to model these observations of the torus; 
most of that w ork pictures the torus as made up 
of clouds (e.g.. iNenkova et al.ll2008a), while soin e 
consider a continuous structure (e.g., Krolikl2007 ) . 



although most of the models are static. In the 
quasi- clumpy models of lDullemond fc van Bemmel 
(l2005l) . "clumps" represented by axisymmetric 
rings of material are used to understand the dif- 
ferences between smooth and clumpy models. As- 
sumptions of their model include the idea that all 
clumps have equal optical depth and a Gaussian 
density profile. They found that several of their 
smooth and clumpy models were able to suppress 
the 10 /Ltm silicate feature in emission. Thus they 
note that although their results corroborate the 
idea that clumpy tori can account for this obser- 
vational feature, it does not unequivocally confirm 
the idea of a clumpy torus, as families of continu- 
ous winds models can explain it as well. 



The formalism of iNenkova et al.l ( 2008a ) takes 



into account the recent interferometric and spec- 
troscopic results. Their model, composed of 
5—15 dusty clouds along each radial equatorial 
ray, is successful at explaining several puzzling as- 
pects of AGN infrared observations. Specifically, 
dust on the bright side of one of their optically 
thick clouds is much hotter than on the dark side, 
and dust on the dark side of a cloud nearer the 
source of the illuminating continuum can be as 
hot as dust on the bright side of one further out. 
The clumpy model furthermore reproduces the be- 
havior of the 10 /im silicate feature, namely the 
lack of very deep absorption features. 



Clu mpy-torus models are clearly useful; iSchartmann et al 
( 20081 ) found with their clumpy tori models that 
the existence of the 10 fim silicate feature, 
whether in emission or absorption, depends on 
the size, optical depth, and distribution of clouds 
closest to the nucleus. However, the family of 
cloud (or clumpy) torus models described above 
typically neglects the dynamics of the clouds; 
a notable exception to this has been work on 
accreting clouds (see e.g 



Vollmer et al.l 2004: 



Beckert fc Dusch]||2004 ISchartmann et al.|[2oTll) . 

In addition, some researchers have postulated 
that these clouds are en trained in a wind (e.g.. 



Elitzur fc Shlosmanll2006[ ). 



1.1. Dynamical Models of the Dusty "Torus" 

The above-mentioned static models of the torus 
do not address the question of why the torus has 
the structure it does, how it lofts gas and dust 
up to large heights above the accretion disk, or 
why its covering fraction may depend on lumi- 
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nosity. To address these questions, a robust, dy- 
namical model must be developed that can ac- 
curately reproduce the features that we observe 
in AGNs. Such dynamical models have been ap- 
proached in a number of ways. O i ie sce nario, first 
proposed by iKrolik &: BegelmanI (Il988l) . suggests 
that the torus consists of a number of optically 
thick clumps which orbit the central engine and 
collide regularly with other clumps. 



In a different dynamical paradigm. iDopita et al 
()l998l) . aiming to describe the IRAS colors of 



Seyfert galaxies, suggested that the torus is a 
large-scale accretion flow of some continuous 
medium towards the nucleus. This free-falling "en- 
velope" of material rotates slowly and circularizes 
once it reaches the centrifugal radius. Since their 
lower limits on the accretion rates are well above 
that which would support the Eddington lumi- 
nosity of the central engine, they determined that 
much of the infalling material wo uld flow away 
from the accretion disk in a wind. iHopk ins et al.l 
(j201l[ ) have recently investigated this kind of dy- 
namical model in numerical simulations, examin- 
ing larger-scale gas flows as a possible model for 
the torus. 

Finally, another possibility is that the torus is a 
dusty wind flowing from the accretion disk. Such 
an idea is supported by the growing evidence for 
gaseous outflows in many types of AGNs. Both 
radio-loud and radio-quiet AGNs often exhibit 
blueshifted absorption lines that can be broad 
and/or narrow; these features are evidence of out- 
flowing material. Approximately 15% of radio- 
quiet quasars display strong, blueshifted absorp- 
tion features at UV resonance transitions with ve- 



locities as high as .1c (e.g.. iReichard et al.ll2003 : 
Gibson et al. I I2OO9I) . There has also been obser- 
vational evidence that suggests the mass outflow 
rate in AGN s is nearly equal to the mass inflow 
rate (see e.g., Crenshaw et al. 20031 Chartas et al 
20031) 



Konigl fc Kartid (1994) first proposed that the 



torus can be explained as a disk-driven hydro- 
magnetic wind. They were motivated to con- 
sider such winds due to their apparent success 
at explaining particular radiative characteristics 
of young stellar objects (YSOs), which are simi- 
lar in some respects to those of AGNs. For ex- 
ample, they note that both types of objects of- 
ten exhibit flat IR spectra, strong Ca III triplet 



lines, and broad Na D emission - observational 
findings which have been interpreted in terms of 
a strong central UV source surrounded by a disk- 
like, dusty mass distribution. The hydro- n iagnetic 
wind model was later expanded upon bv lEverett 
(I2OO5I) to include a more realistic treatment of ra- 
diative acceleration in magnetic winds. 

In addition, models where infrared radia- 
tion pushes a continuous wind vertically off of 
the accretion d isk bey ond the dust sublima- 
tion radius (e.g., 'Krolik'2007 t IShi fc Krolikll2008t 
Dorodnitsvn et al. 2011ai ;bl) have also been ex- 
amined. This vertical, infrared radiation pressure 
may indeed be present, but the radiation force due 
to the central, UV-bright AGN continuum has no- 
tably more integrated power to drive dust grains 
than the IR continuum; whereas, for instance. 



Dorodnitsvn et al.l (|2011b[ ) find a radiation pres- 
sure on dust from infrared photons can increase 
the effective Eddington ratio by a factor of approx- 
imately 10-30, our calculations show that, when 
including the entire AGN continuum, that factor 
can become 200 (this is somewhat less, but com- 
parable to the factor o f ^ 800 enhancement f ound 
in the early models o f lKonigl fc Kartielll994l) . 

This paper will focus on the dynamical model 
of the torus as a dusty wind, which is launched 
by both magnetohydrodynamical (MHD) forces 
as well as radiation pressure due to the accretion 
disk continuum. The dusty wind generated in this 
manner can cover a large fraction of the sky as seen 
from the central black hole. This model has the 
benefit of being more self-consistent and includ- 
ing the important physics of motion around the 
black hole and radiative acceleration. Ultimately, 
we wish to use this model to understand how the 
physical properties of dusty winds in AGNs cor- 
relate with their observable spectral signatures in 
the IR. 

2. The Structure of the Torus and the 
Wind Model 



Our model of the torus consists of a dusty 
wind driven by both hydro-magnetic forces and 
radiative acceleration. The model and its corre- 
sponding code has been e xpanded upon from its 
ori ginal form (deta iled bv lKonigl fc Kartig|l994l ) 



by lEverettj (|2005l ). where a comprehensive ac- 



count of the model's components and key equa- 
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tio ns are describ ed. In this model, we advance 
Everettl (|2005l ) by adding the continuum opac- 



on 



ity of ISM dust grains, as specified by the ISM 
dust mode l in Cloudy (version 06.02.09b, last de- 
scr ibed by Ferland et al. 1998t for th e dust model. 



see 



Mathis et al.lll977l: Ivan Hoof et a l. 2004) 



A s in previous magneto-cent rifugal models (see, 
e.g., iBlandford fc Pavnd Il982l hereafter BP82), 
the model assumes a parsec-scale magnetic field, 
approximately vertical to the accretion disk. (It 
is important to point out that the origin of such 
a field is not clear, although simulations of jet 
launching from AGN accretion disks seem to favor 
the advection of magnetic flux from large scales; 
this is a topic of great interest in current research. 



see e.g.. lHawlevll2011l ) With such a field, gas and 
dust particles orbiting in the accretion disk are 
subject to the centrifugal force, pointing outwards 
along the disk, as well as the gravitational force 
directed inwards towards the central black hole. 
However, a charged particle, tied to the magnetic 
field, will be flung outward along the magnetic 
field line if the centrifugal acceleration overtakes 
the gravitational acceleration; this can occur if the 
angle of the magnetic field line to the vertical is 
> 30° (BP82). In this sense, the "foundation" of 
our model is a straightforward a pplication of the 
self- si milar model of BP82 and Konigl fc Kartjd 
(jl994[ ). Such self-similar models allow relatively 
quick calculations of the wind geometry and dy- 
namics by intrinsically assuming that the shape of 
the magnetic field lines (which change due to radi- 
ation pressure) scales as radius, such that the wind 
geometry at large distances has a similar shape to 
the wind geometry at small radii, only scaled up 
to a larger size. This self-similarity is a significant 
assumption, but it allows for the calculation of 
both the radial and vertical momentum equations, 
and in this case, allows for the simple addition of 
the radiative acceleration by effectively decreasing 
the g ravitational potential (see the Appendix in 
Ever ett 2005). 

As part of the assumption of self-similarity, the 
model requires that the density at the surface of 
the accretion disk scales as a function of radius, so 
that po oc r^^ and Bq oc Tq '^^^^^/'^ ^ where the sub- 
script '0' denotes values at the disk surface. Our 
particular implementation of the wind model al- 
lows for any radial scaling parameter, 5, but in all 
of the calculations presented here, we set 6 = 3/2 



(this is equivalent to assuming a stationary accre- 
tion disk flow, where the mass-loss does not change 
across different decades in disk radius) . Two other 
key parameters that are set in this model are: A, 
the ratio of total specific angular momentum in 
the wind to that in the disk, and k, the dimen- 
sionless ratio of mass flux to magnetic flux; as in 
the 'standard' model of BP82, we set these param- 
eters to A = 30 and k — 0.03. These parameters 
help set the magnetic field strength in the wind at 
the surface of the disk; as opposed to jet solutions 
of BP82 which required B ~ 10^ to 10* G, the 
strongest field in our models was of order 0.03 G 
at the base of the wind. 

Once the wind is launched by centrifugal ac- 
celeration, it is also subject to radiative acceler- 
ation from the photons emitted from the accre- 
tion disk. The radiation pressure on the dust 
from the central source is very strong - even for 
cases where L/LEdd = 0.1, it is approximately 
10 times the force required to unbind dust from 
the gravitational potential ( Everett et al. 2009[ ). 
(Though line scattering is included in the calcu- 
lation, the continuum source dominates the line- 
radiation pressure in the dusty medium.) Dust 
grains absorb radially (anisotropically) streaming 
photons originating in the accretion disk, and then 
de-excite isotropically. The force, felt by the dust 
particles due to conservation of momentum, feeds 
back on the wind structure by bending the mag- 
netic wind radially away from the source. The ra- 
diative force therefore works in conjunction with 
the magneto-centrifugal forces to accelerate the 
wind flow off the accretion disk, and modifies the 
structure of the outflow. 

The assumed dust composition in the model is a 
standard interstellar distribution, with a mixture 
of silicates and graphites and a continuous distri- 
bution of grain sizes. These dust grains will be 
destroyed at temperatures greater than their sub- 
limation temperature, Tg^b 1500 K. The wind's 
innermost radius for dust driving is therefore set 
approximately to the dust sublimation radius, out- 
side of which the silicate and graphite dust can 
survive. The sublimation radius is dependent on 
a number of factors, including the bolometric lu- 
minosity and shape of the incident continuum. 

Schartmann et al. I (|2008h note during their 
investigations of the dusty torus as a clumpy 
medium that it is important to have a variety of 
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sublimation radii for different sizes of dust grains 
in order to accurately model the IR SED from the 
dust. Their solution was to utilize a model with 
three different grain species, each with five dif- 
ferent grain sizes. Such a multi-component grain 
distribution would be important, but given the 
complexity of our models and the large amount 
of time required for computing each SED in even 
the simplified single-grain case, we will use a sin- 
gle grain type and sublimation radius as our first 
approximation. 

2.1. The Model Code 

The model code works in two separate steps. 
First, the structure of the wind is determined semi- 
analytically, starting with a pure MHD wind and 
then adding in radiation pressure from dust and 
atomic lines. Second, a Monte Carlo simulation 
uses the information about the structure of the 
wind to generate a spectral energy distribution 
from 0.1 — 2000 /im. As we are interested primarily 
in the more readily observable IR regime, we plot 
these results from 2 — 100 fiui. Furthermore, at 
wavelengths < 2 /zm, we cannot be certain of the 
accuracy of our SED calculations, as the accuracy 
is affected by the parameters chosen for the Monte 
Carlo simulation (see § 13. 1|) . 

In more detail, the program starts by calcu- 
lating the pure m agneto-centrifu gal wind dynam- 
ics (see Figure [U lEverett 2005[) without includ- 
ing the effects of radiation pressure. Parame- 
ters such as the black hole mass, column den- 
sity at the base of the wind and starting loca- 
tion of the wind are specified here. In this phase, 
an estimation of the density and velocity struc- 
ture of the wind as a function of height is calcu- 
lated by a relatively simple semi-analytic model of 
the magneto-centrifugal acceleration of g as aloni 
a stre amline (this is d escribed in detail in lEveret 
20051) . Next, Cloudy (|Ferland et al.lll998l ) is used 
to calculate the photo-ionization of the gas and de- 
termine opacities throughout the structure of the 
wind, given the central continuum. With all of this 
information, the radiation pressure on the wind is 
calcul ated. This is exactly as described in Everett 
POOSI ). except for the addition of the opacity of 
dust, which is supplied via Cloudy's "ISM" dust 
model. This information is passed back into the 
first step, where the MHD wind structure is modi- 
fied accordingly. This process is iterated until con- 



vergence, which typically takes seven iterations. 

Shown in Figure [2] are the gas streamlines for 
our fiducial model (with parameters as specified 
in § HI) , plotted in the z direction along the axis 
of symmetry as a function of the radial (r) direc- 
tion. The streamlines are shown for each iteration 
of the MHD radiative wind code. The radiation 
pressure acts to bend the streamlines radially out- 
ward, away from the central source, with each it- 
eration, until the code converges. 

The second stage of our calculatio n uses the 
Monte Carlo simulation program MC3D (jWolfll2003h 
to predict spectral energy distributions. The in- 
put for the simulation is the output from the MHD 
radiative wind code, which includes the struc- 
ture of the wind and the SED of the central con- 
tinuum source (we use the composite SED from 
Richards et al.ll2006L modified to remove IR emis- 
sion; see Fig. [3] and S 13.31 for more details). This 
input is passed to MC3D: this is a general purpose, 
3D Monte Carlo code which calculates heating and 
emission when given a central radiative source sur- 
rounded by a scattering and absorbing medium. 
MC3D first determines the temperature of the dust 
grains due to heating from the continuum and ree- 
mission from other dust grains. In a second step, 
MCSD's ray-tracer is then used to create the SEDs 
from 0.1 - 2000 ^m, which we plot from 2 - 100 
/im. 

3. Methodology 



Initial tests for basic functionality, stability and 
consistency of the code are documented in detail in 
Everett (|2005h . For this work, the code was com- 
piled on a number of different machines and it was 
verified that identical results were produced on 
each machine. With confidence that the code per- 
forms as expected, we then investigated ways to 
minimize computation time while maintaining an 
acceptable level of accuracy. The bulk of the com- 
puting time is spent on generating the IR SED, 
and thus we analyzed the MC3D input parameters 
to determine both whether we could make our sim- 
ulations more accurate without a significant in- 
crease in time, and whether the process could be 
sped up significantly without a loss of accuracy. 
These tests are reviewed in § 13.11 

After these basic considerations, we set out 
to explore the physical parameter space that the 
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model allows. AGNs exhibit a large range of ob- 
served parameters, with varying luminosities, cen- 
tral continuum shapes, central black hole masses, 
etc., and the way these parameters interact in pro- 
ducing an IR SED is not necessarily simple or obvi- 
ous. Salient parameters were identified and tested 
against one another in order to determine their ef- 
fects on the dusty wind, and therefore on the shape 
and power of the IR SED. An in-depth description 
of each parameter examined can be found in § 13.21 
to §11 

3.1. MC3D Testing Parameters 

MC3D is a versatile code, with a range of cus- 
tomizable user-set parameters for modelling dust- 
temperature distributions and SEDs. We tested 
a number of these in order to determine an ap- 
propriate set of base parameters which maintain 
a reasonable degree of accuracy while minimizing 
the amount of time needed to generate the models. 

MC3D has a number of model geometry options 
to choose from; because our model is axisymmet- 
ric, we utilized the fully two-dimensional model 
with a radial and vertical dependence for the den- 
sity distribution (see Figure S]). Our default base 
values for testing the effects of varying these pa- 
rameters were as follows: a maximum model ra- 
dius. Rout, of 20 parsecs, with 30 sub-divisions 
(grid points) in the radial direction and 101 sub- 
divisions in the 9 direction, a half-opening angle, 
ip, of 75 degrees, and 100 photon packets per wave- 
length. We will discuss in more detail, below, our 
tests of these parameter settings; for all of the pa- 
rameter changes discussed. Figure [5] displays the 
differences in the SED and Table [1] details the ef- 
fects on the emission at 2, 3, 6, 10, and 50 /im as 
compared to the default parameters. 

For example, decreasing the number of photon 
packets per wavelength from the default of 100 to 
10 (Wolf .200 3) had little effect on the SED, but 
decreased the time required to compute the tem- 
perature distribution by almost a factor of five. 
The variation that is seen between SEDs generated 
with 10 or 100 photon packets per wavelength is 
a minor loss of accuracy in the emission from the 
hottest dust; the greatest difference is at 2 /im, 
where the model with 10 photon packets shows 
1.61 times the emission as the model with 100 
photon packets. As this loss is not relevant for our 
purposes - the largest discrepancy is seen at wave- 



lengths < 3 /im where the accretion disk and host 
galaxy typically contribute significantly to the ob- 
served emission - and the increase in computation 
efficiency is great, we opted to use 10 photon pack- 
ets per wavelength for all of our simulations. 

Increasing the resolution of the grid by choosing 
201 6* subdivisions increased the computation time 
for the SED by several hours, with little benefit. 
Again, the largest difference is seen at 2 fim, where 
the test model with 201 6* subdivisions showed 1.11 
times the emission compared to the default of 101 
9 subdivisions. 

Decreasing Rout to 10 parsecs had little effect 
on the time for generating the temperature dis- 
tribution, but increased the time for generating 
SEDs by more than a factor of 3. The difference 
in the SEDs is slight. The greatest difference is 
seen at long wavelengths; at 50 /tm, the test model 
with Rout = 10 pc displays 0.88 times the emission 
compared to the default model with Rout = 20 pc. 
We retain the value of 20 pc for Rout ■ 

Decreasing the half-opening angle from ip — 75° 
to ^ — 60° also served to increase, rather than 
decrease, the computation time. The SED showed 
a very slight overall change in normalization, with 
the model with -0 = 60° showing a minor decrease 
in emission. 

To summarize, for the rest of our models, we 
take the MC3D parameters as: 10 photon packets 
per wavelength, 101 9 sub-divisions. Rout = 20 pc, 
and V = 75°. 

A full list of these parameters, the range of val- 
ues explored, and the corresponding times to gen- 
erate a temperature distribution and SED can be 
found in Table [5] Most of the computations were 
done on a workstation with a dual-core CPU, each 
at 3.16 GHz and a workstation with a quad-core 
CPU, each at 2.83 GHz. 

3.2. Exploring the Parameter Space 

Turning to more physical parameters, we have 
investigated variations in the final IR SED due 
to inclination angle of the observer (i), the base 
column density of the wind {Nh,o), and changes 
to the central continuum shape, luminosity, and 
black hole mass. All parameter changes are spec- 
ified before the wind structure is calculated, with 
the exception of the inclination angle, which is 
only specified when MC3D is called. We chose to run 
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many of our models at an inclination of i = 60°, 
as choosing inclination angles smaller than that 
tended to increase the computation time. Table |3] 
details these parameters along with corresponding 
relevant information, such as the terminal veloc- 
ity of the wind and the mass outflow rate, as cal- 
culated by the MHD wind model, as well as the 
kinetic luminosity of the outflows. 

The relative power and shape of the SED were 
only slightly affected by a change in the inclina- 
tion angle, with larger inclination angles showing a 
slightly lower normalization in the overall IR SED 
(see Figures [6] and [71). This indicates that the IR 
emission is consistent with being optically thin at 
most wavelengths of interest. 

The integrated luminosity varied significantly 
with changes in the column density, with higher 
columns producing a higher overall infrared lumi- 
nosity. We chose to maintain a fixed column den- 
sity Nh,o = 10^^ cm~^ for the remaining models, 
as this was the only column which could match 
the order of magnitude of luminosity comparable 
to a typical quasar (L w 10"*^ erg s~^). Such a 
model generates mass outflow rates on order of 
IMq yr-i (see Table EJ. 

3.3. The Central Continuum 

In order to investigate what sort of radiative 
output we would expect to see from the dusty 
wind, we specify a central input continuum that 
represents an approximation of the emission from 
the accretion disk around the black hole. We used 
a composite created from broadband photometry 
of 259 optically luminous Sloan Dig ital Sky Survey 
(SPSS ) quasars, as described by [Richards et al 



(|2006l ). as a starting point for this central contin- 
uum (see Figure |8] for a plot of this continuum, as 
well as the input continua modified from it). The 
optical and UV portions of the spectrum are as- 
sumed to be a fair representation of the emission 
from the accretion disk. However, the observed 
"IR hump" (from 1 fim - 100 fim) in this com- 
posite spectrum is due to the dust emission that 
we aim to model. Therefore, the IR hump was 
removed and replaced with a simple power law ex- 
trapolated from the optical continuum. Though 
the observed IR hump from 1 - 100 /im is quite 
prominent, common features such as the 10 and 
18 /im silicate emission features are not apparent, 
due to a lack of spectral resolution in the broad- 



band photometry and the effect of smearing out 
such features with composite averaging. Never- 
theless, such a composite SED serves as an excel- 
lent basis for comparing to our IR SEDs models, 
as a check to ensure that the power and general 
shape can be matched. We will be testing whether 
the wind model can account for the empirical IR 
hump. 

3.3.1. Changes to the Central Continuum: aox 

The spectral index aox describes the amount 
of energy emitted in the X-ray relative to the 
amount emitted in the UV/optical, and can be 
considered a measure of the "X-ray brightness" 
of a source. The pararn eter is defined as follows 
(jTananbaum et al.l[l979l ): 



-0.384 log 



L^{2 keV) 



L^(2500 A) 



(1) 



Quasars have been found with a range of aox 
values, typically varying from 1.2 to 1.8. There is 
an anti-correlation between aox and UV luminos- 
ity, as more luminous AGNs show less X-ray emis- 
sion r elative to emission in the optical wav elengths 



e.g.. lJust et al.ll2007t ISteffen et al.ll2006[ ). In or 



der to test this range of aox values, we selected 
an "average" value of 1.6 (appropriate for the UV 
luminosity of our fiducial model), and then modi- 
fied the input continuum (see Figure [Sj to reflect 
extreme values of aox = 1-1 (much higher than av- 
erage output in the X-ray), and aox = 2.1 (much 
lower than average output in the X-ray). The X- 
ray modification begins at 7 x 10^^ Hz and extends 
to 10^^ Hz. The discontinuity at = 7 x 10^^ Hz 
is not physical, but serves to isolate the effect of 
the X-ray continuum from that of the far-UV. 

These new continua are passed as parameters 
to the wind and radiative acceleration code. The 
results of these changes are described in § ID 

3.3.2. Changes to the Central Continuum: L/ LEdd 

Initial tests of the wind model were performed 
with an Eddington ratio of L/ LEdd ~ 0.01. Such 
an Eddington ratio is characteristic of the Seyfert 
regime, but may not accurately describe quasar 
luminosities. In order to test the model at the 
higher luminosities more appropriate to quasars, 
we incrementally increased the Eddington ratio to 
L/LEdd = 0.1 by increasing the total luminosity 
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and keeping Mbh fixed. Furthermore, we tested 
the effects of changing black hole mass, which has 
the effect of changing the Eddington luminosity. 



smce: 



^Edd 



1.25 X fO' 



38 



M 



erg s 



(2) 



We tested masses ranging from Mbh = 10^ — 
10^ Mq. We note that although the part of the 
code that generates the wind structure was able 
to accommodate such high luminosities with rel- 
ative ease, the time for generating the final SED 
increased by a factor of three or more at higher 
luminosities. 

4. Results 

As a basis of comparison for many of our re- 
sults, we adopt a fiducial model with the following 
parameters: Nh^ = 10^^ cm~^, L/Lsdd = 0.1, 
Mbh = 10* Mq,\ = 60°, and aox = 1.6. The Ed- 
dington luminosity for a black hole of mass AIbh 
= 10* Mq is 1.25 X lO'*^ erg s^^, giving a luminos- 
ity at 5100 A of L5100 = 3.7 x lO^^ erg s'^ Hz"! 
for an Eddington ratio of L/Lsdd = 0.1. Figure |8] 
shows this fiducial model, plotted with the input 
continuum, and the original Richards et al. ( 2006[ ) 
composite SED. The sum of the predicted IR SED 
and the input SED is also displayed. The compos- 
ite SED is normalized to the input SED at Lsiooj 
and the output SED has been integrated over Att 
steradians to match the units of the observations. 

With the base column density of Nh,o = 10^^ 
cm^^, the predictions ap proximately ma t ch the 



power expected from the iRichards et al. ( 2006h 
composite SED. The observed composite SED has 
~ 82% as much power radiated in the IR (from 
2-100 ^m) as in the optical/UV (from 2 ^m- 1000 
A) When normalized at isiooi our model with 
the smaller column density of Nh.o — 10^^ cm~^, 
L/LEdd = 0.01, Mbh = 10* M©, i = 60°, and 
Q-ox = 1-6 has only 45% as much power radiated 
in the IR as in the composite optical/UV, which 
does not match the expected output. Our model 
with a higher column of Nh,o = 10^^ cm~^ (and all 
other parameters as before) radiates 97% as much 
power in the IR as in the composite optical/UV, 
which is more than what is expected. Our fiducial 
model, with a higher column of Nh,q = 10^^ cm^^ 
and a higher Eddington ratio of L/LEdd = 0.1, has 



~ 70% as much power radiated in the IR as in the 
composite optical/UV, which is closer to the ex- 
pected 82%. 

While the model SEDs can account for the gross 
power and shape of the composite SED with the 
parameter values in the fiducial model, there are 
specific features that do not match. Specifically, 
the composite has more power at both shorter 
(2-8 /im) and longer (> 30 /im) wavelengths, and 
the fiducial model SEDs show more prominent 
and sharply peaked silicate features at both 10 
and 18 /im. On this second point, note that the 
composite by construction smears out spectral 
structure such as the silicate emission bumps as 
it is made using broad-band photometry; many 
mid-IR spectra of luminous quasars, however, 
show prominent silicate emission features (see 



ow p: 

e.g.. ISiebenmorgen et al.l ,2005: Hao et al.l 1200 



J 



Netzer et alTl2007l : iMason et al.ll2009t iDeo et al 



20111 ). We also note that at long wavelengths, we 
do not expect much of an AGN contribution from 
the dusty wind. Empirically, far-infrared emis- 
sion of AGNs is typically weak unle ss th e host 
is actively star-fo rming (Schweitzer et al. 2006t 



Netzer et al.ll2007^ . In fact, the composite SED 



contains little data beyond 24 /im, with that part 
of the SED constructed by "gap repair," replacing 
the missing data with information extrapolated 
from lElvis et al.l ( 19941 ) who used IRAS photom- 
etry of low redshift quasars at these wavelengths. 
We explore these issues in more detail in § [5l 

4.1. Mass outflows, terminal velocities, 
and kinetic luminosities 

Table |3] contains information on the mass out- 
flow rates (dM/dt), terminal velocities (wterm), 
and kinetic luminosities (Lkm) of each of the mod- 
els. We see mass outflow rates ranging from 
0.96 — 4.I9M0 yr~^, terminal velocities ranging 
from 1904 — 7992 km sec"""^, and kinetic luminosi- 
ties ranging from 1.13 x 10^^ — 7.71 x 10*'^ erg 
s~^. These values are strongly affected by the 
black hole mass of the model: the largest mass 
outflows, terminal velocities, and kinetic luminosi- 
ties are produced by the models with the highest 
black hole masses, even at lower Eddington ratios. 
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4.2. Results of varying inclination angle 

Here, we examine the result of varying the in- 
clination angle from i = 90° (edge-on) to z = 30° 
(nearly face-on). Inclination angles smaller than 
i = 30° were not considered because of the great 
increase in time required to compute the SEDs. As 
the first phase of our parameter exploration, these 
models were run at luminosities L/LEdd — 0.01 
and 0.03, lower than that of our fiducial model. 

Figure [6] shows the result of varying the in- 
clination angle for a model with N^.o = 10^''' 
cm-2, L/LEdd = 0.01, Mbh = 10** Mq, and 
ctox = 1-6. This model is different from our fidu- 
cial model in two ways: it has both a lower col- 
umn density {Nh.o — lO^'* cm~^, as compared to 
Nh,o = lO^'^ cm^^), and a lower Eddington ratio 
{L/LEdd — 0.01, as compared to L/LEdd — 0.1). 
TablelUdetails the amount of power emitted in the 
IR as compared to the optical/UV for each model, 
as well as the ratios of emission at 3, 6, 10, and 
50 /im, relative to the i = 90° model. Varying 
the inclination angle for these parameters has al- 
most no effect on the shape or normalization of 
the SEDs. The lack of variation indicates that IR 
emission must be optically thin at the wavelengths 
of interest. 

However, when we move to a higher base col- 
umn density, variations in the SED start to be- 
come apparent. Figure [7| shows the effects of incli- 
nation angle for a model with Nh,o — 10^^ cm~^, 
L/LEdd = 0.01, Mbh = lO^M©, and = 1.6. 
This model is different from our fiducial model as 
it has a lower Eddington ratio {L/LEdd = 0.01, 
compared to L/LEdd = 0.1). Variations in the 
SED are much more apparent than they were at 
the lower column density of Nh,o = 10^*^ cm~^. 
Table [5] details the amount of power emitted in 
the IR as compared to the optical/UV for each 
model, as well as the ratios of emission at 3, 6, 10, 
and 50 /im, relative to the i = 90° model. 

At wavelengths > 5 fim, the effect of decreasing 
the inclination angle is primarily a change in nor- 
malization, with smaller inclination angles (more 
face-on lines of sight) displaying more power. 
At wavelengths < 5/im, the effect is more pro- 
nounced, and the SEDs display a change in shape 
as well as normalization, with larger inclinations 
(more edge-on lines of sight) showing a marked de- 
crease in the amount of emission from the hottest 



dust. At 3/im, the i = 30° SED displays 2.4 
times more emission than the i = 90° SED, and 
1.48 times more emission at 50 /xm. The results 
of varying inclination angle are similar when look- 
ing at the model with a slightly higher Eddington 
ratio, L/LEdd — 0.03, which is still smaller than 
that of our fiducial model. 

The measurable variation at higher column 
densities indicates that the wind is becoming more 
optically thick — especially at the shortest wave- 
lengths — as the column increases. 

4.3. Results of varying aox 

Next, we examine the effects of changing the 
amount of energy emitted in the X-ray relative to 
the amount emitted in the UV/optical, defined by 
parameter aox- Figure |9] displays the results of 
varying aox as compared to our fiducial model, 
with iV^^,o = 10^5 cm-2, L/LEdd = 0.1, Mbh = 
10^ M0, and i ~ 60°. For reference, a Attox value 
of 1.0 corresponds to a factor of 400 difference 
in X-ray luminosity at 2 keV for a set value of 
L - . The depressed X-ray emission of the X- 

2500 A ^ / 

ray faint model crudely mimics the effect of strong 
X-ray absorption seen in some broad a bsorption 
line quasars fe.g.. lGallagher et al.|[2006h . The X- 
ray bright model with aox = 1.1 displays the high- 
est overall luminosity at wavelengths > 4 /im, with 
emission that decreases at wavelengths < 4 /im, 
relative to the other models. Our fiducial model, 
with aox = 1-6, shows somewhat lower emission 
than both the X-ray bright and X-ray faint mod- 
els. 

Table [S] shows the amount of power emitted in 
the IR as compared to the optical/UV for each 
model, as well as the ratios of emission at 3, 6, 
10, and 50 /im, r elative to the ap^ = 1 .6 fiducial 
model. When the iRichards et all (200a) compos- 



ite is normalized to the ^5100 value of the input 
continuum, the fiducial model displays 69.7% as 
much power in the I R as compared to the opti- 
cal/UV power in the Richards et al. ( 20061 ) com- 
posite SED. The aox = 1-1 and aox = 2.1 models 
display 102% and 84.9% as much power in the IR 
as compared to the composite optical/UV, respec- 
tively. 

The X-ray faint model with aox = 2.1 displays 
more emission at shorter wavelengths (< 4/im) 
than either of the other models, and at wave- 
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lengths > 4 /im shows a change in normahzation 
compared to the other models. At wavelengths 
> 50 /J,m, the SEDs for all the models appear to 
converge. Compared to the fiducial model with 
dox = 1.6, the greatest difference is seen for the 
dox = 1-1 model which shows 1.55 times more 
emission at 6 /im. 

These results are perhaps counterintuitive, and 
reveal the complex role of the SED in both accel- 
erating and heating the dust. The X-ray bright 
model has significantly more power available to 
illuminate the dusty medium, and therefore it 
makes sense that the IR SED is more luminous 
relative to the other two models. However, the rel- 
ative lack of emission at the shortest wavelengths 
in the X-ray bright model indicates that there is 
a smaller volume of dust at the highest tempera- 
tures, a natural consequence of a wind that is more 
radially "combed out" (and hence has a narrower 
wedge of material near the sublimation tempera- 
ture) . 

In addition, the dust sublimation radius, i?subj 
is ~ 13% larger for the aox=l-l models relative 
to the aox=l-6 models (see TablelHl). We consider 
the 2% difference in sublimation radii between the 
aox=l-6 and 2.1 runs to be the same, within the 
precision of the simulation. 

These changes highlight the complex nature 
of how enhanced X-ray emission impacts winds. 
First, while X-rays contribute to heating, their 
efficiency is quite low compared to UV photons: 
at 300 eV (near the Oox jumps), the dust opac- 
ity (scattering plus absorption) is about an order 
of magnitude b elow the UV opacity (see Fig. 9 
m lDrainel[2003h . Therefore, it requires a strong 
increase in the X-ray flux to yield, e.g., a factor 
of two increase in the near-IR flux. A second ef- 
fect is the increased momentum transfer per pho- 
ton of X-rays vs. UV light: X-rays are both more 
likely to scatter (versus being absorbed) and carry 
more momentum. Supporting this interpretation, 
the X-ray- faint model with aox = 2.1 shows the 
strongest short (< 4 /im) wavelength emission. 
The UV power provides a "floor" to i?sub; the 
smaller X-ray power in the aox = 2.1 model does 
not change the sublimation radius significantly rel- 
ative to the aox = 1.6 model, while increasing the 
relative X-ray power by setting aox = 1-1 prefer- 
entially pushes out i?sub and accelerates the wind 
more rapidly. Both effects reduce the volume of 



dust heated to the highest temperatures, and thus 
the relative power coming out at the shortest IR 
wavelengths. 

Figure [TU] displays the results of varying ctox for 



a model with ( 



10 



25 



Edd 



0.01, 



Msff = 5 X 10^ Mq , and i = 60° . This model has a 
lower Eddington ratio than our fiducial model, as 
well as a higher black hole mass. The X-ray bright 
model with ctox = 1.1 displays the highest overall 
luminosity. The model with ctox = 1.6 shows a 
change in normalization, with less emission than 
the aox = 1.1 model, and little effect on the SED 
shape. Table [7] shows the amount of power emit- 
ted in the IR as compared to the optical/UV for 
each model, as well as the ratio of emission at 3, 
6, 10, and 50 /xm, relative to the aox = 1-6 model. 
When the Richards et al. ( 2006h composite is nor- 
malized to the isioo value of the input continuum, 
the model with aox = 1.6 displays about 80% as 
much power in the I R as compared to the opti- 
cal/UV power in the Richards et al.l ([2006) com- 
posite SED. The aox = 1-1 and aox = 2.1 models 
display ^ 120% and ~ 90% as much power in the 
IR as compared to the composite optical/UV, re- 
spectively. 

The X-ray faint model with aox — 2.1 also 
shows a change in normalization compared to the 
X-ray bright model, with less overall emission, and 
is very similar to the model with aox = 1-6, though 
it displays a slight increase in emission at wave- 
lengths < 8 /im. Compared to the model with 
Ciox = 1-6, the aox — 1-1 model shows 1.75 times 
more emission at 3 /im and 1.27 times more emis- 
sion at 50 /im, and the aox = 2.1 model shows only 
1.20 times more emission at 3 /im and displays the 
same amount of emission at 50 /im. As in the pre- 
vious set of models, the impact of changing aox 
from 1.6 to 1.1 is to push out i?sub by 18%; the 
difference between the two X-ray fainter models of 
2% is not considered significant. 

The difference in the wind response to varying 
aox models when the values of Mbh and L/LEdd 
are also changed illustrates the interplay between 
these parameters, and the challenges in isolating 
the effects of the SED. An increase in Mbh of a 
factor of 5 coupled with a decrease in L/LEdd of a 
factor of 10 means that the luminosity is a factor 
of 2 fainter than in the previous model described 
in this section. A lower luminosity coupled with a 
larger black hole mass generates a wind with more 
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vertical streamlines, because the net radial com- 
ponent of the ejection (centrifugal plus radiative 
minus gravitational) force will be reduced. For ex- 
ample, in the L/L^dd — 0.1 X-ray bright model, 
the wind streamlines at large radii have an angle 
from the disk of 63° vs. 74° for the lower luminos- 
ity, X-ray bright model. The more vertical wind 
structure of the lower luminosity models keeps 
the dusty wind closer to the illuminating source 
for longer, but apparently decreases the observed 
emission from the hottest dust. We consider it 
likely that this is an artifact of the observer's in- 
clination angle (see Figure [7|) ; the more vertical 
(because of lower radial acceleration) wind will be 
more optically thick at shorter wavelengths. The 
view of the hottest dust is therefore obscured for 
an observer viewing angle of 60°. 

4.4. Result of varying the black hole mass 

We examine the effects of changing the mass 
of the central black hole, with values from Mbh 
= 10* - 10^ Mq. Altering Mbh also changes the 
Eddington luminosity (see Equation [2]). For Mbh 
= 10* M0, the Eddington luminosity is Lsdd = 
1.25 X 10^^ erg s^^, which gives a luminosity at 
5100 A of L5100 = 3.7 X 1028 erg s'^ Hz"! for 
L/LEdd = 0.01. The luminosities scale linearly 
with the black hole mass. In these tests, we keep 
L/LEdd = 0.01 fixed. 

Figure [TT] shows the effect of changing the black 
hole mass, which changes the luminosity for a 
fixed L/LEdd- In this situation, higher black hole 
masses produce overall more luminous SEDs. For 
both AfB// = 5 X 10* Mq and Mbh = 10^ Mq, 
the shape and normalization of the SED is affected 
from 2—40 /im as compared to the Mbh = 10* Mq 
model. Table [8] shows the amount of power emit- 
ted in the IR as compared to the optical/UV for 
each model, as well as the ratios of emission at 
3, 6, 10, and 50 ^m, relative to Mbh = 10* Mq. 
The higher Mbh models show a more peaked SED 
shape from 3 — 30 /^m as compared to the Mbh = 
10* M0 model, as well as the overall increase in 
luminosity. Compared to the model with Mbh = 
10* Mq, the model with Mbh = 5 x 10* Mq dis- 
plays ^ 60 times the emission at 10 /xm and ^ 80 
times the emission at 50 /zm, whereas the model 
with Mbh = 10^ Mq displays ^ 130 times the 
emission at 10 /xm, and ^ 40 times the emission at 
50 /im. Given that the emission at these compar- 



ison wavelengths is not just scaling linearly with 
the luminosity (which is increasing by a factor of 5 
and 10 as Mbh is increased by the same factors), 
the change in black hole mass clearly has a direct 
influence on the IR power. 

4.5. Result of varying L/LEdd 

We examine the effects of changing the Edding- 
ton ratio, L/LEdd, which, for a given black hole 
mass, has the effect of changing the bolometric lu- 
minosity. For Mbh = 10* Mq, an L/LEdd ratio 
of 0.1 gives a luminosity at 5100 A of L5100 = 
3.7 X 10^9 erg s^^ Hz-i; L/LEdd = 0.07 gives 
L5100 = 2.6 X 1029 erg s-i Hz-i; L/LEdd = 0.03 
gives isioo = 1-1 X 10^9 erg s^^ Hz"^; and an 
L/LEdd = 0.01 gives igmo = 3.7 x 10^* erg s'^ 
Hz-i. 

Figure [12] displays the effects of varying L/LEdd 
(changing the luminosity). As is expected, higher 
L/LEdd ratios display a greater luminosity as com- 
pared to lower L/LEdd ratios. The L/LEdd = 0.1 
and 0.07 models are very similar in shape, differ- 
ing only by a normalization factor. The L/LEdd = 
0.03 and 0.01 models show a marked change in 
SED shape as compared to the higher luminosity 
models, in addition to displaying a lower overall 
luminosity. Table [9] shows the amount of power 
emitted in the IR as compared to the optical/UV 
for each model, as well as the ratios of emission 
at 3, 6, 10, and 50 /im, relative to L/LEdd — 0.1. 
As compared to the fiducial model with L/LEdd — 
0.1, the L/LEdd — 0.07 model displays around 0.8 
times the emission at each wavelength examined, 
whereas the L/LEdd = 0.03 model displays ~ 0.4 
times the emission at 3 /im and 0.5 times the emis- 
sion at 50 /im. 

4.6. Distinguishing between Mbh effects 
and L/LEdd effects 

As changing both Mbh and L/LEdd effectively 
results in a change in luminosity, one might ask 
whether the changes to these parameters are re- 
dundant. To isolate these effects, we plot two 
models with the same input luminosity, £5100 = 
3.7 x 10^9 erg s'^ Hz"! (see Figure The 
two plots are very similar, but not identical: the 
model with a higher black hole mass and lower Ed- 
dington ratio displays a higher normalization from 
2 — 40 /im. These results suggest that the prop- 
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erties are not redundant, though running further 
models for comparison would be of great benefit. 

5. Conclusions and Discussion 

Wc have investigated the parameter space of 
our model of the torus as a dusty magneto- 
centrifugal wind with radiative acceleration. We 
have examined a number of interesting and im- 
portant parameters and illustrated how they af- 
fect the IR SED of an AGN. We have verified 
that with input SEDs of bolometric luminosities 
L K, IQ'^'^ — 10^^ erg s^"'^ and a base column den- 
sity of 10^^ cm~^, the code produced reasonable 
IR SEDs with approximately the right shape and 
luminosi ty {L ^ 10^'^ ~ I Q"^^ ^.' "S s~^) as expected 



from the IRichards et al.l (|2006l ) composite of opti- 
cally luminous SDSS quasars. This is a promising 
result for a relatively simple model given that we 
have not attempted any fitting. A benefit to our 
model is that we are able to see directly the ef- 
fects of various physical parameters on the final 
IR SED, which ultimately will allow us to under- 
stand the physical properties of the torus itself. 
By determining which physical parameters have 
an observable effect on the IR SEDs, and nar- 
rowing down these parameters to ascertain which 
ones allow us to generally reproduce the power 
expected, we have established a reasonable start- 
ing point from which we can expand and further 
refine our model. 

We summarize our results as follows: 

1. Varying the inclination angle has little ef- 
fect at base column densities < 10^^ cm~^, 
but becomes more pronounced at base col- 
umn densities > 10^^ cm~^. This indicates 
that the emission is becoming slightly more 
optically thick as the column increases. The 
most salient difference (see Fig. 8) is a deficit 
at the shortest wavelengths for the larger 
inclination angles, indicating that emission 
from the hottest dust, at the smallest radii, 
is being blocked. 

2. The parameter aox , characterizing the amount 
of energy emitted in the X-ray relative 
to the amount emitted in the optical/UV, 
has slight but measurable effects on the 
IR SEDs. The X-ray bright model, with 
c^ox = 1-1, displays the highest overall lu- 



minosity, with emission that decreases at 
wavelengths < 4 /xm. Both the X-ray bright 
(ctox = 1-1) and X-ray faint {uox — 2.1) 
models display marginally more emission 
than the fiducial model with aox = 1-6. 

3. The short wavelength (< 5 ^m) emission ap- 
pears to be sensitive to the relative power of 
X-rays relative to UV in the SED, offering 
promise for using the properties of the IR 
SED to constrain the strength and shape of 
the high energy continuum. This illustrates 
the utility of models where the IR-emitting 
medium can respond dynamically to the in- 
put SED; in static models the incident spec- 
trum is quickly thermalized and information 
on its shape is lost. 

4. Higher black hole masses produce more lu- 
minous IR SEDs, for a fixed value of L/LEdd- 
The SEDs for M^jj > 5 x lO^Mg are more 
peaked from ~ 3 — 30 ^m as compared to the 
SEDs for Msff < 5 X 10^ Mq. 

5. Higher Eddington ratios produce more lumi- 
nous IR SEDs. The SEDs for L/LEdd > 0.03 
are more peaked from ^3 — 30 /im as com- 
pared to SEDs for L/LEdd < 0.03. 

6. Although changes to both L/LEdd and Mbh 
are effectively changes to the overall input 
luminosity, these effects may not be redun- 
dant, and appear to be slightly distinguish- 
able from each other. The model with a 
higher black hole mass and lower Eddington 
ratio displays a higher normalization from 
2-40 ^m. 

At present, although our model is producing 
approximately the right amount of power, it does 
not adequately reproduce the amount of expected 
emission in the hot dust, at wavelengths < 8 /im. 
This may be due in part to inclination angle ef- 
fects. As discussed in § 14. 2[ although inclination 
angle has little effect at lower column densities, 
it begins to show an effect at column densities of 
^H,o = 10^^ cm~^, with smaller inclination an- 
gles displaying more power, particularly at wave- 
lengths < 5 fim. It is possible that this effect is 
amplified not only by an increase in column den- 
sity, but by an increase in luminosity as well. We 
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chose an inclination angle oi i — 60° for the ma- 
jority of our models in order to increase our com- 
putational efficiency, which allowed us to run more 
models in a shorter period of time. However, such 
an inclination angle is not likely to represent the 
average for Type 1 objects; a more characteris- 
tic choice would have been between i — 45° and 
i = 30°. The choice of z = 60° means that we may 
be obscuring the innermost parts of the wind, and 
therefore missing the hottest dust, in many of our 
models. 

Another possibility is that our assumption of an 
ISM dust composition is not entirely correct. We 
assume a single silicate dust sublimation radius, 
which neglects the hot graphite grains that can 
exist at smaller radii, and thus we are missing that 
emission. Further investigations into the details of 
our dust composition and sublimation radius are 
therefore required. Changing grain composition 
can also affect the profile of the silicate features. 

The above simulations with MC3D required long 
run-times. The fiducial model takes ~3-4 days 
on a 4-core desktop machine to construct an SED 
for a single viewing/inclination angle. Further, 
smaller sublimation radii and smaller inclination 
angles require much longer run-times. For exam- 
ple, our fiducial model with a modified i?sub = 
1.5675 X 10^^ cm requires '--^2-3 months on an 8- 
core desktop machine. Smaller sublimation radii 
translate to higher grid resolutions in the grid 
setup procedure for MC3D. To explore the sensitiv- 
ity of our results to the effects of a smaller sublima- 
tion radius, we found it necessary to port MC3D to 
run on the many-core symmetric multi-processing 
(SMP) compute nodes of SHARCNET0. The fidu- 
cial model now takes 19 hours when using 32 
nodes; the highest resolution run mentioned above 
takes slightly longer than a week. 

We modified our fiducial model, setting the 
inner radius of the dusty wind to be 1.5675 x 
10^* cm, about half that of the original fiducial 
model, keeping all other parameters fixed. A max- 
imum temperature of 1615 K was reached within 
the dusty wind as modeled by MC3D, as opposed 
to approximately 800-1200 K in previous varia- 
tions of the fiducial model. This model repro- 
duces the blackbody-like shape of the near-IR 



emission between 1-8 ^m (see Figure [T4|) that is 
seen in mid-I R spect ra of luminou s quasars (e.g., 
Netzer et al.l i,2007: D eo et al.ll201l[ ). The SED cor- 



responding to our fiducial run (red curve) and the 
high-resolution run (blue curve), bracket the av- 
erage near-IR emission in the composite SED of 
Richards et al. (2006). The composite mean SED 
is an average over several important parameters 
such as inclination, luminosity and Mbh, and for 
a more accurate comparison, one would have to 
average over the similar range of models. That 
full study is beyond the scope of this paper, but 
we find that, in our predicted SEDs, the near-IR 
hot-dust luminosity decreases as we increase the 
inclination angle, giving approximately the cor- 
rect trend. These results show that changes in the 
dust sublimation radius approximate the expected 
presence of hot graphite grains, and corresponding 
grid resolution (an effect of the choice of grid con- 
struction process) are important effects to consider 
when comparing model SEDs to observations. 

We have explored a wider parameter space of 
this model based ori ginally on t he hydro-magnetic 
wind model of iKonigl &: Kartie (1994). Our re- 
sults are roughly consistent with theirs, with our 
models displaying similar 10 /xm micron emission 
to th eir model (labeled as 'A2' in Konigl fc Kartid 
19941 ). with Mbh = 10^ Mq, a bolometric lumi- 
nosity Lboi = lO''^ erg s""'^, and a mass outflow 
rate of 1.0 yr~^. We do not, however, see the 
silicate feature in absorption, as they do with their 
model 'A3', with Mbh = lO'^ M©, a bolometric lu- 
minosity Lfioi = 10^"* erg s^^, and a mass outflow 
rate of 3.0 M© yr^^, even in our models with ex- 
tremely high mass outflow rates ( > 4.0 Mf7^ yr~^) 
and overall higher mass. Nenko va et al. I (l2008bl) 



^ SHARCNET is a consortium of Canadian academic institu- 
tions who share a network of high performance computers. 



also see the silicate feature in absorption with 
many of their clumpy tori models. This may in 
part mean that a higher resolution is required to 
resolve the wind, although models with a signifi- 
cantly increased resolution would also require an 
increase in computa tion time (see above di scus- 
sion). Furthermore, Konigl &: Kartie ( 1994f ) add 
in radiative acceleration only approximately; it is 
possible that with more realistic radiative acceler- 
ation added in, the wind accelerates off the disk 
quickly enough that the density drops very quickly 
with height, and only a very small region at the 
disk surface remains where we might observe sili- 
cate absorption. 
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Curiously, we find that t he shape of our 



SEDs, in general, are similar to iSchartmann et al 



( 20081 ) 's clumpy, face-on models, and are rather 
dissimilar to their continuous models (see the up 



per l eft panel in Figure 10 of ISchartmann et al 



20081 ) . We assume th at they take the same par am- 
eters as described in ISchartmann et alJ (j2005[ ): a 
black hole mass of Mbh — 6.6 x lO^'M©, and 
an Ldisk/LEdd value of 0.06, parameters which 
are descriptive of the lower luminosity Seyfert 
regime of AGNs. Although our models are geared 
more towards describing luminous quasars, their 
model is most comparable to our model with 
Mbh = 1 X 10* M© and L/LEdd = 0.07. Again, 
however, we do not see the silicate feature in ab- 
sorption as they do for edge-on models. 

5.1. Future Work 

Future work on this model will entail refining 
our parameters, such as testing smaller inclina- 
tion angles to determine whether we can more ac- 
curately reproduce the emission from the hottest 
dust. Along these lines, we aim to decrease the 
computation time further, to enable us to resolve 
the inner graphite-dominated region that likely 
gives rise to the near-lR emission. Then, we will 
be able to explore the importance of multiple dust 
components, at the necessary higher resolution. 
Further next steps will involve detailed compar- 
isons of the predicted SEDs with a broad sample 
of IR observations, for example from the Spitzer 
Space Telescope. Such in-depth comparisons will 
allow us to further constrain the salient parame- 
ters that we have isolated in this paper. 

Currently, we are investigating the effects of in- 
creasing the "shielding gas" column density, Ngg. 
The high X-ray flux from the AGN is prone to 
over-ionizing the gas in the wind. A thick layer of 
hot, highly ionized "shielding gas," located inte- 
rior to the wind-launching dust sublimation ra- 
dius, is required in order to prevent this over- 
ionization. Such a construct was first hypothe- 
sized in the context of continuous q uasar winds as 



the "h itchhiking gas" describe d bv iMurrav et al 



1995 ). an d later s t udied by IProga fc Kallman 
(120041 ) and lEverettl (|2005l ). Subsequent empiri- 



cal evidence for the shielding gas was found in 
the signatures of high column density X-ray ab- 
sorption in broad ab sorption line quasars (e.g. 
GaUagher et al.ll2006[) . Modifications to the MC3D 



code are required in order to deal with absorbed 
continua passing through different column densiti- 
ties of shielding gas to determine the effect on the 
dust sublimation radius and output IR SED. 

Finally, another step towards perhaps more re- 
alistic models would be to simulate dusty winds 
of a filamentary nature. Such a model may bet- 
ter represent the structures seen i n interferomet 
ric ob servations of AGNS (see, e.g 
20071) . 
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Table 1 

Variations in testing parameters for MC3D 



Parameter Change 2 /im 


Ratio 3 /xm Ratio 6 fim Ratio 


10 /xm Ratio 


50 fim Ratio 


Default 


1.00 


1.00 1.00 




1.00 


1.00 




10 photons/A 1.61 


1.14 0.99 




0.99 


0.98 




201 divisions 1.11 


1.04 1.00 




0.99 


1.00 




Rout = 10 


pc 1.04 


1.04 1.02 




0.99 


0.88 




tjj = 60° 


0.98 


0.96 0.96 




0.96 


0.96 




IN 1 bj. — . 


All models are compared to the model with default base parameters of 100 photon 


packets per 


wavelength, 101 


6 sub-divisions, Rout = 20 pc, and ip 


= 75°. 










Table 2 














Testing Parameters 








Model 


Photon Packets Model Radius # of 9 






i 4 


tsED^ 


Name 


per A 


(pc) Divisions 


n 


n (cm-2) 


(hrs) 


(hrs) 


n4testSl 


100 


20 101 


75 


60 10^4 


38 


59 


n4testS2 


10 


20 101 


75 


60 1024 


8 


58 


n4testS3 


100 


20 201 


75 


60 1024 


39 


64 


n4testS4 


100 


10 101 


75 


60 1024 


38 


187 


n4testS5 


100 


20 101 


60 


60 1024 


49 


77 


n5 tests 1 


10 


20 101 


75 


30 1025 


4 


31 


n5testS2 


10 


20 101 


75 


45 1025 


4 


37 


n5testS3 


10 


20 101 


75 


60 1025 


4 


34 


n5testS4 


10 


20 101 


75 


75 1025 


4 


31 


nStestSS 


10 


20 101 


75 


90 1025 


4 


44 



^Half-opening angle 
2 Inclination angle 
^Base column density 

4Time to generate temperature distribution 
5 Time to generate SED 
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Table 3 
Parameter Space 



Model 


Mbh' 


L/ LEd<P 




dM/df^ 
(Mo yr-i) 


'^term 


r 6 


Name 


(xlO^ M0) 






(km sec~^) 




n5elx0s6 




0.01 


1.6 


0.99 


1904 


1.13 


n5elxas6 


J 


0.01 


1.1 


0.96 


2344 


1.67 


n5elxbs6 




0.01 


2.1 


0.99 


1959 


1.19 


n5e2x0s6 




0.03 


1.6 


1.00 


2755 


2.40 


n5e2xas6 




0.03 


1.1 


1.05 


2884 


2.76 


n5e2xbs6 




0.03 


2.1 


1.00 


2754 


2.39 


n5e3x0s6 




0.07 


1.6 


1.14 


2877 


2.97 


n5e4x0s6 




0.1 


1.6 


1.23 


2899 


3.25 


n5e4xas6 




0.1 


1.1 


1.32 


3162 


4.14 


n5e4xbs6 




0.1 


2.1 


1.22 


2907 


3.24 


n5elinls6 


5 


0.01 


1.6 


2.39 


5994 


27.1 


elnilxas6 


5 


0.01 


1.1 


2.51 


5958 


28.1 


elmlxbs6 


5 


0.01 


2.1 


2.38 


6000 


27.0 


n5e4inls6 


5 


0.1 


1.6 


3.80 


4999 


29.9 


e4nilxas6 


5 


0.1 


1.1 


4.19 


5268 


36.6 


e4mlxbs6 


5 


0.1 


2.1 


3.78 


5006 


29.9 


n5elm2s6 


10 


0.01 


1.6 


3.82 


7992 


76.9 


elm2xas6 


10 


0.01 


1.1 


4.07 


7758 


77.1 


elni2xbs6 


10 


0.01 


2.1 


3.82 


7991 


76.8 



Note. — All models have a base column density of Nnfi = 10^^ cm~^. 
^ Black hole mass 
^Eddington ratio 

^Amount of power in optical relative to X-ray 
*Mass outflow rate 
^Terminal velocity of outflow 
^Kinetic luminosity of outflow 
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Table 4 

Variations in i for Nha = 10^"* cm~^ 



Parameter Change i^iRji^optlvv 3 /im Ratio 6 /im Ratio 10 /im Ratio 50 /ini Ratio 



i = 90° 44.3% 1.00 1.00 1.00 1.00 

i = 75° 46.0% 1.03 1.04 1.04 1.11 

i = 60° 45.5% 0.96 1.04 1.02 1.12 

i = 45° 47.9% 1.02 1.07 1.09 1.14 

i = 30° 49.8% 1.07 1.09 1.15 1.15 



Note. — All niodels are compared to the model with i — 90°, N^^i^ = 10^^ cm^^, L/LEdd — 
0.01, Mbh = 10* Mq, and ttox = 1.6. The ratio Lm/ L^pi/uy is the amount of power emitted in 
the IR from 2 — 100 /im by a certain model as compared to the amount of power em itted in the 
optical/UV (from 2 /im — 1000 A) by the composite continuum bv lRichards et al. (2006), when both 
are normalized at isioo- 



Table 5 

Variations in i for Nha = 10^^ cm~^ 



Parameter Change Lm/L^pi/uy 3 /im Ratio 6 /im Ratio 10 /im Ratio 50 /im Ratio 



i = 90° 70.9% 1.00 1.00 1.00 1.00 

i = 75° 85.3% 1.59 1.14 1.14 1.34 

i = 60° 97.5% 1.91 1.30 1.31 1.42 

i = 45° 106.6% 2.13 1.39 1.47 1.46 

i = 30° 114.7% 2.40 1.48 1.65 1.48 



Note. — All models are compared to the model with i — 90°, Nh,o = 10 cm^ , L/LEdd — 
0.01, Mbh = 10* M©, and aox — 1-6. The ratio Lm/ L^pifuy is the amount of power emitted in 
the IR from 2 — 100 /im by a certain model as compared to the amount of power em itted in the 
optical/UV (from 2 /im — 1000 A) by the composite continuum bv lRichards et al" (2006), when both 
are normalized at L^iqq. 
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Table 6 

Variations in aox for L/LEdd = 0.1 and Mbh = 10^ M© 



Parameter Change 


LirJ Lopt/uv 


Rsnh (pc) 


3 ^Tn Ratio 


6 /im Ratio 


10 /im Ratio 


50 /im Ratio 


aox = 1-1 


101.9% 


1.78 


1.49 


1.55 


1.54 


1.14 


aox = 1.6 


69.7% 


1.57 


1.00 


1.00 


1.00 


1.00 


aox = 2.1 


84.9% 


1.54 


1.56 


1.34 


1.18 


1.09 



Note. — The ratio Lm/ Lopt/uv is the amount of power emitted in the IR from 2 — 100 /im by a certain model 
as compared to the amou nt of power emitted in the optical/UV (from 2 //m— 1000 A) by the composite continuum 
by [Richards et al~ ( 2006 ). when both are normahzed at Lsioo- In columns 4 through 7, all models are compared 



to the model with a„ 



1.6, 



60°, N. 



H,0 



10 



25 



cm 



'-,L/L 



Edd 



0.1, and Mbh = 10** M©, 



Table 7 

Variations in aox for L/LEdd = 0.01, and Mbh = 5 x 10® Mq 



Parameter Change 


Lir/ Lopt/uv 


-Rsub (pc) 


3 fiia Ratio 


6 fim Ratio 


10 /im Ratio 


50 /im Ratio 


aox = 1.1 


123.8% 


0.58 


1.75 


1.63 


1.50 


1.27 


aox = 1.6 


83.3% 


0.49 


1.00 


1.00 


1.00 


1.00 


aox = 2.1 


88.0% 


0.50 


1.20 


1.13 


1.04 


1.00 



Note. — The ratio Lir/ Lopt/uv is the amount of power emitted in the IR from 2 — 100 /im by a certain model 
as compared to the amount of power emitted in the optical/UV (from 2 /im— 1000 A) by the composite continuum 
by lRichards et al" (2006), when both are normalized at L5100. For columns 4 through 7, all models are compared 
to the model with aox = 1.6, i = 60°, iV^.o = 10^^ cm'^, L/LEdd = 0.01, and Mbh = 5 x 10® M©. 



Table 8 
Variations in Mbh 



Parameter Change 


Lir/ Lopt/uv 


3 /im Ratio 


6 /im Ratio 


10 /im Ratio 


50 /im Ratio 


Mbh = 10® Mq 


97.5% 


1.00 


1.00 


1.00 


1.00 


Mbh = 5 X 10® M© 


83.3% 


3.23 


4.71 


6.13 


1.94 


Mbh = lO'^ M© 


82.6% 


7.46 


10.35 


12.55 


2.96 



Note. — All models are compared to the model with Mbh = 10® M©, aox = 1.6, i = 60°, Nh.q — 
10^^ cm~^, and L/Lsdd = 0.01. The ratio Lir/ Lopt/uv is the amount of power emitted in the IR from 
2 — 100 /im by a certain model as compared to t he amount of power em itted in the optical/UV (from 
2 /im — 1000 A) by the composite continuum bv lRichards et al" ( 20061 ). when both are normalized at 
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Table 9 
Variations in L/LEdd 



Parameter Change 


LibJ Lgpi/jjy 


3 /im Ratio 


6 /im Ratio 


10 /xm Ratio 


50 /im Ratio 


L/LEdd = 0.1 


69.7% 


1.00 


1.00 


1.00 


1.00 


L/LEdd = 0.07 


79.6% 


0.80 


0.79 


0.78 


0.84 


L/LEdd - 0.03 


87.9% 


0.36 


0.35 


0.35 


0.51 


L/LEdd = 0.01 


97.5% 


0.16 


0.12 


0.096 


0.36 



Note. — All models are compared to the model with L/Lsdd = 0.1, Mbh = 10® M©, aox = 1.6, 
i = 60°, and Nh,o = 10^^ cm^^. The ratio Lm/ Lopt/uv is the amount of power emitted in 
the IR from 2 — 100 /im by a certain model as compared to the amount of power em itted in the 
optical/UV (from 2 /im — 1000 A) by the composite continuum bv lRichards et al~ (2006), when both 
are normalized at L5100. 



20 




Central AccretioiiDisk 
Black Hole 



Fig. 1. — Schematic diagram of the structure 
of the radiat ively accelerated magneto-centrifugal 
wind model (jEver ctt 2005). The sub parsec-scale 
accretion disk is the source of the continuum illu- 
minating the wind. The wind has two distinct 
components. The leftmost, thinner component 
is the "shielding gas," and is modeled as a pure 
magneto-centrifugal wind. The wider streamline 
on the right is the dusty wind which is propelled by 
both magneto-centrifugal and radiative accelera- 
tion. The black lines within the shield and stream- 
line are examples of the zones in which Cloudy de- 
termines opacities and simulates photo-ionization. 
In our calculations, the shielding gas is set to have 
an effectively negligible column density so that the 
wind is dominated by the radiatively accelerated 
component. 



Gos Streamlines: Fiducial Mode 




Fig. 2. — Gas streamlines for the fiducial model 
with Nh^o = 10^5 cm-2, L/LEdd = 0.1, Mbh - 
10® Mq, aox — 1-6, and i = 60°. Streamlines are 
plotted in the z direction (along the axis of sym- 
metry) in units of parsecs, as a function of the r 
(radial) direction, also in units of parsecs. The 
streamlines are shown for each iteration of the 
MHD radiative wind code. At each subsequent 
iteration (until convergence), the wind becomes 
more radial as radiative acceleration is taken into 
account. 
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[Hz]) 



Fig. 3. — The input continuum showing the 
range of a. ox values used. Luminosity is shown 
in units of log vLu, and frequency is in units of 
log Hz. The solid l i ne sho ws the original SED 



from Richards et al 



( 200d ). which includes the 
IR hump. The dotted, dashed and dashed-dotted 
lines display the modified continuum, with the IR 
hump replaced by a power law, and with an aox 
of 1.1, 1.6, and 2.1 respectively. The discontinu- 
ity at j/ = 7 x 10^'^ Hz (or 0.29 keV), where the 
X-ray continuum modification begins, is not phys- 
ical, but serves to isolate the effect of the X-ray 
continuum from that of the far-UV. 



Fig. 4. — MC3D model geometry. Listed in 



Wolj (|2003l ) as geometry 2D(b), it is a fully two- 



dimensional model, with vertical and radial depen- 
dences for the density distributions. The model is 
axisymmetric about the z axis. Shown are the 
MC3D grid divisions in the Q and r directions. 
The half-opening angle, t/;, is measured from the 
horizontal. The model radius is defined by Rout- 
The inclination angle, z, is defined from the nor- 
mal, such that i — 90° is perpendicular to the axis 
of symmetry. 
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N„^(,= 10"* cr^"", L/Le,, = 0.01, a„=1.6, i = 60" 




2 10 100 

A [Mm] 



Fig. 5. — Testing the effects of variation in the 
MC3D parameters. These IR SEDs are generated 
by MC3D with an inchnation angle of z = 60°, a 
base column density of A'^^.o — 10^^ cm~^, an Ed- 
dington ratio of L/LEdd — 0.01, a black hole mass 
of Mbh = 10® Mq, and an aox value of 1.6. Lu- 
minosity is in units of vL^, and wavelengths are 
in /im. Here, "default" refers to the set of base 
testing parameters in § 13. li namely a model ra- 
dius. Rout, of 20 parsecs with 30 sub-divisions in 
the radial direction, 101 sub-divisions in the 9 di- 
rection, a half opening angle, tjj, of 75 degrees, and 
100 photon packets per wavelength. Each subse- 
quent SED varies one of these parameters, leaving 
the rest fixed at the default values. The only sig- 
nificant decrease in computation time arose from 
the change from 100 to 10 photons per wavelength, 
which is also the only result that deviated slightly 
from the SED generated by the base parameters. 
The largest discrepancy is seen at wavelengths 
< 3 /Ltm, where the accretion disk and host galaxy 
will also contribute to the observed emission. 




2 10 100 



Fig. 6. — Result of varying the observer's inclina- 
tion angle, i, for a model with Nh^ = 10^* cm~^, 
L/Lsdd = 0.01, Mbh = IO^Mq, and a^, = 1.6. 
Luminosity is shown in units of vL^^ and wave- 
length in units of /im. The SEDs display remark- 
ably little variation with changes in the observer's 
inclination angle. 
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Fig. 7. — Result of varying the observer's inclina- 
tion angle, i, for a model with Nh^ = 10^^ cm~^, 
L/LEdd = 0.01, Mbh = IO^Mq, and = 1.6. 
Luminosity is shown in units of I'L^, and wave- 
length in units of /um. At wavelengths > 5/im, 
the effect of decreasing the inclination angle is pri- 
marily a change in normalization, with smaller in- 
clination angles displaying more power. At wave- 
lengths < 5 /im, the effect is more pronounced, 
and the SEDs display a change in shape as well as 
normalization, with larger inclinations showing a 
steeper decrease in the amount of observable emis- 
sion from the hottest dust. 
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Fig. 8. — Fiducial model with iV^f_o = 10^''' cm ^, 
L/LEdd = 0.1, Mbh = IO^Mq, Uox = 1.6, and 
i = 60°. Luminosity is shown in units of log uLi^, 
and wavelength in units of log Hz. The output 
SED for the fiducial model is denoted by the solid 
line. Also shown is the input SED (dotted line), 
the sum of the output and input SEDs (dashed-dot 
line), and the Richards composite SED (dashed 
line). 
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N„„=10^= cm"", L/L„, = 0.1, i = 60, M„„=10° M,, 




N„,,= 10^^ cm"^ L/Lej, = 0.01, i = 60', Mg„ = 5x10° M^, 







Fig. 9. — Result of varying Uox for a model with 
Nh,o = 1025 cm-2, L/LEdd = 0.1, Mbh = 
10^ Mq, and i = 60°. Luminosity is shown in 
units of i>L,^, and wavelength is in units of /im. 
The model with Uox = 1-1 (solid line) displays the 
highest overall luminosity, with emission that de- 
creases at wavelengths < 4 /im. The model with 
= 1-6 (dotted line) displays a slightly lower 
emission than both the X-ray brighter and X-ray 
fainter models. The model with ~ 2.1 (dashed 
line) has a similar SED shape to the other models, 
but displays more emission at shorter wavelengths 
(< 4/im) than cither of the other models. The 
shape and emission for wavelengths > 50 /xm is 
very similar for all models. 



Fig. 10. — Result of varying aox for a model 
with Nh,o = 10^5 cm-2, L/LEdd = 0.01, Mbh 
= 5 X 10* Mq, and i = 60°. Luminosity is shown 
in units of vLi,, and wavelength is in units of /um. 
The X-ray bright model with aox = 1-1 (solid line) 
displays the highest overall luminosity. The model 
with = 1-6 (dotted line) shows primarily a 
change in normalization, with little effect on the 
SED shape. The X-ray faint model with aox = 2.1 
(dashed line) also shows a change in normalization 
compared to the X-ray bright model, and is very 
similar to the model with aox = 1-6, though it dis- 
plays a slight increase in emission at wavelengths 
< 8/um. 
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Nh,o=10 cn \ L/Lejj = 0.01, i = 60, a„=1. 




N„„=10^= cm"^ i = 60 , a„,= 1.6, Mb„=10° Mj^ 




Fig. 11. — Result of varying Mbh for a model with 
Nh,o = 10^^ cm-2, L/Lsdrf = 0.01, a^^ = 1.6, 
and i = 60°. Luminosity is shown in units of vL^, 
and wavelength is in units of /um. The highest 
black hole mass, Mbh = 10® (dashed line), 
shows an increase in himinosity at all wavelengths 
compared to the model with Mbh = 5 x 10^ Mq 
(dotted line) and Mbh = 10^ Mq (solid line). 
Here, higher black hole masses produce higher 
overall luminosities. The higher Mbh models 
show a more peaked SED shape from 3 — 30 /xm as 
compared to the Mbh = 10^ Mq model, as well 
as the overall increase in luminosity. 



Fig. 12. — Result of varying L/LEdd for a model 
with Nh,o = 10^^ cm-2, Mbh = 10^ M©, aox = 
1.6, and i = 60°. Luminosity is shown in imits 
of vLi,, and wavelength is in units of /^m. The 
difference between L/LEdd = 0.1 (solid line) and 
L/LEdd = 0.07 (dotted line) is a normalization 
factor, with little change in the SED shape. More 
pronounced changes to the SED shape, in addi- 
tion to a change in normalization, can be seen for 
L/LEdd = 0.03 (dashed line) and L/Lsdd = 0.01 
(dashed-dot line), with the latter displaying rel- 
atively more emission than the former at wave- 
lengths < 6/im. 
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Fig. 13. — Plot of models with the same input 
luminosity, L5100 = 3.7 x 10^^ erg s~^ Hz^^, for 
different Eddington ratios and black hole masses. 
The solid line shows our fiducial model, with 
Mbh = 10** M0 and L/LEdd = 0.1. The dot- 
ted line shows a model with Mbh = 10^ and 
L/LEdd = 0.01. The two curves are not iden- 
tical, with the model with Mbh = 10^ Mq and 
L/LEdd — 0.01 showing relatively more power 
shortward of the 10/im silicate feature. 
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log u [Hz] 



Fig. 14. — The fiducial model with a smaller dust- 
sublimation radius shows an SED (blue solid line) 
peaking in the near-IR. As the dust-sublimation 
radius is increased the SED power decreases pro- 
gressively in the near-IR and the peak emission 
shifts to longer wavelengths. 



REFERENCES 

Antonucci, R. 1993, ARA&A, 31, 473 

Antomicci, R. R. J., & Miller, J. S. 1985, ApJ, 
297, 621 

Barvainis, R. 1990, ApJ, 353, 419 

Beckert, T., & Duschl, W. J. 2004, A&A, 426, 445 

Blandford, R. D., & Payne, D. G. 1982, MNRAS, 
199, 883 

Chartas, G., Brandt, W. N., & Gallagher, S. C. 
2003, ApJ, 595, 85 

Crenshaw, D. M., Kraemer, S. B., & George, 1. M. 
2003, ARA&A, 41, 117 

Deo, R. R, Richards, G. T., Nikutta, R., et al. 
2011, ApJ, 729, 108 

Dopita, M. A., Hcislcr, C., Lumsden, S., & Bailey, 
J. 1998, ApJ, 498, 570 

Dorodnitsyn, A., Bisnovatyi-Kogan, G. S., & Kail- 
man, T. 2011a, ApJ, 741, 29 

Dorodnitsyn, A., Kallman, T., & Bisnovatyi- 
Kogan, G. S. 2011b, ArXiv e-prints 

Draine, B. T. 2003, ARA&A, 41, 241 

Dullemond, C. P., & van Bemmel, I. M. 2005, 
A&A, 436, 47 

Elitzur, M., & Shlosman, I. 2006, ApJ, 648, LlOl 

Elvis, M., Wilkes, B. J., McDowell, J. C., et al. 
1994, ApJS, 95, 1 

Everett, J. E. 2005, ApJ, 631, 689 

Everett, J. E., Gallagher, S. C., & Keating, S. K. 
2009, in American Institute of Physics Con- 
ference Proceedings, Vol. 1201, The Monster's 
Fiery Breath; Feedback in Galaxies, Groups, 
and Clusters, ed. S. Heinz & E. Wilcots, 56-59 

Ferland, G. J., Korista, K. T., Verner, D. A., et al. 

1998, PASP, 110, 761 

Fritz, J., Franceschini, A., & Hatziminaoglou, E. 
2006, MNRAS, 366, 767 

Gallagher, S. C, Brandt, W. N., Chartas, G., 
et al. 2006, ApJ, 644, 709 



Gibson, R. R., Jiang, L., Brandt, W. N., et al. 
2009, ApJ, 692, 758 

Grimes, J. A., Rawlings, S., & Willott, C. J. 2004, 
MNRAS, 349, 503 

Hao, L., Weedman, D. W., Spoon, H. W. W., et al. 
2007, ApJ, 655, L77 

Hao, L., Strauss, M. A., Fan, X., et al. 2005a, AJ, 
129, 1795 

Hao, L., Spoon, H. W. W., Sloan, G. C, et al. 
2005b, ApJ, 625, L75 

Hawley, J. F. 2011, in lAU Symposium, Vol. 275, 
lAU Symposium, ed. G. E. Romero, R. A. Sun- 
yaev, & T. Belloni, 50-58 

Hill, G. J., Goodrich, R. W., & Depoy, D. L. 1996, 
ApJ, 462, 163 

Honig, S. F., & Kishimoto, M. 2010, A&A, 523, 
A27-h 

Hopkins, P. F., Hayward, C. C Narayanan, D., 
& Hernquist, L. 2011, MNRAS, 2115 

Jaffe, W., Meisenheimer, K., Rottgering, H. J. A., 
et al. 2004, Nature, 429, 47 

Just, D. W., Brandt, W. N., Shemmer, O., et al. 
2007, ApJ, 665, 1004 

Konigl, A., & Kartje, J. F. 1994, ApJ, 434, 446 

Krolik, J. H. 2007, ApJ, 661, 52 

Krolik, J. H., & Begelman, M. C. 1988, ApJ, 329, 
702 

Lawrence, A. 1991, MNRAS, 252, 586 

Lawrence, A., & Elvis, M. 2010, ApJ, 714, 561 

Mason, R. E., Levenson, N. A., Shi, Y., et al. 2009, 
ApJ, 693, L136 

Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 
1977, ApJ, 217, 425 

Murray, N., Chiang, J., Grossman, S. A., & Voit, 
G. M. 1995, ApJ, 451, 498 

Nenkova, M., Sirocky, M. M., Ivezic, Z., & Elitzur, 
M. 2008a, ApJ, 685, 147 



28 



Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, 
Z., & Elitzur, M. 2008b, ApJ, 685, 160 

Netzer, H., Lutz, D., Schweitzer, M., et al. 2007, 
ApJ, 666, 806 

Pier, E. A., & Krolik, J. H. 1992, ApJ, 401, 99 

Proga, D., & Kallman, T. R. 2004, ApJ, 616, 688 

Reichard, T. A., Richards, G. T., Hall, P. B., et al. 
2003, AJ, 126, 2594 

Richards, G. T., Lacy, M., Storrie-Lombardi, 
L. J., et al. 2006, ApJS, 166, 470 

Sanders, D. B., Phinney, E. S., Neugebauer, G., 
Soifer, B. T., & Matthews, K. 1989, ApJ, 347, 

29 

Schartmann, M., Krause, M., & Burkert, A. 2011, 
MNRAS, 415, 741 

Schartmann, M., Meisenheimer, K., Camenzind, 
M., Wolf, S., & Henning, T. 2005, A&A, 437, 
861 

Schartmann, M., Meisenheimer, K., Camenzind, 
M., et al. 2008, A&A, 482, 67 



Schweitzer, M., Lutz, D., Sturm, E., et al. 2006, 
ApJ, 649, 79 

Shi, J., & Krolik, J. H. 2008, ApJ, 679, 1018 

Siebenmorgen, R., Haas, M., Kriigel, E., & Schulz, 
B. 2005, A&A, 436, L5 

Simpson, C. 2005, MNRAS, 360, 565 

Simpson, C., & Rawlings, S. 2000, MNRAS, 317, 
1023 

Steffen, A. T., Barger, A. J., Cowie, L. L., 
Mushotzky, R. P., & Yang, Y. 2003, ApJ, 596, 

L23 

Steffen, A. T., Strateva, I., Brandt, W. N., et al. 
2006, AJ, 131, 2826 

Tananbaum, H., Avni, Y., Branduardi, G., et al. 
1979, ApJ, 234, L9 

Tristram, K. R. W., Meisenheimer, K., Jaffe, W., 
et al. 2007, A&A, 474, 837 

Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 



van Hoof, P. A. M., Wcingartncr, J. C., Martin, 
P. G., Volk, K., & Ferland, G. J. 2004, MNRAS, 
350, 1330 

Vollmer, B., Beckert, T., & Duschl, W. J. 2004, 
A&A, 413, 949 

Wolf, S. 2003, Computer Physics Communica- 
tions, 150, 99 



This 2-column preprint was prepared with the AAS lATJ^X 
macros v5.2. 



29 



